function x = cheby_points( a, b, n )

x  = zeros(1,n);
for i = 1:n;

    z =  -cos(( 2 * i - 1 ) * pi / ( 2 * n ));

    x(i) = (z + 1) * (b - a) / 2 + a;
    
end